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Abstract 

A  perturbation  model  for  high-frequency  sound  scattering  from  an 
irregular  elastic  sea  bed  is  considered.  The  sea  bed  is  assumed  ho¬ 
mogeneous  on  the  average  and  two  kinds  of  irregularities  are  assumed 
to  cause  scattering:  roughness  of  the  water-sea  bed  interface  and  vol¬ 
ume  inhomogeneities  of  the  sediment  mass  density  and  the  speeds  of 
compressional  and  shear  waves.  The  first-order  small  perturbation  ap¬ 
proximation  is  used  to  obtain  expressions  for  the  scattering  amplitude 
and  bistatic  scattering  strength.  The  angular  dependence  of  the  scat¬ 
tering  strength  is  calculated  for  sedimentary  rock  and  the  influence 
of  shear  elasticity  is  examined  by  comparison  with  the  case  of  a  fluid 
bottom.  Shear  effects  are  shown  to  be  strong  and  complicated. 

PACS  numbers:  43,30Pc,  43.60.Pt,  43.30.Gv,  43.30Ma 


INTRODUCTION 

The  Born  approximation  (first-order  perturbation  theory)  has  been  applied 
to  sound  scattering  by  rough,  elastic  sea  beds  by  Lapin  [23,  24],  Dacol  and 
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Berman  [4],  and  Kuo  [19].  Kuperman  and  Schmidt  [21]  have  used  the  Born 
approximation  to  obtain  the  coherent  field  in  layered  elastic  sea  beds  with 
rough  interface.  Essen  [5]  has  also  used  the  Born  approximation  to  treat 
rough  sea  bed  types  ranging  from  very  soft  sediment  to  basalt.  This  previous 
work  on  interface  scattering  leaves  open  the  question  of  the  effect  of  shear  on 
scattering  from  within  the  volume  of  the  sea  bed.  Numerous  authors  have 
applied  the  Born  approximation  to  volume  scattering  in  fluid  sediments  [7,  8, 
9,  11,6,  18,  25,  26,  27,  28,  32],  and  these  results  show  that  volume  scattering 
can  be  an  important  mechanism.  Crowther  [3],  Ivakin  [10,  12,  14,  15,  16] 
and  Kuo  [20]  have  applied  the  Born  approximation  to  the  joint  roughness- 
inhomogeneity  problem  for  fluid  sediments.  While  Ivakin  [13]  has  given  a 
formal  Born  approximation  treatment  of  the  volume  scattering  problem  for 
elastic  sea  beds,  there  has  been  no  joint  treatment  of  roughness  and  volume 
scattering  for  elastic  sea  beds.  The  purpose  of  the  present  article  is  to  put 
the  rather  complicated  roughness  and  volume  scattering  theoretical  results 
into  a  common  and  general  format  that  is  straightforward  from  a  numerical 
point  of  view  and  that  lends  itself  to  systematic  study  of  these  two  processes. 
A  subsequent  article  [17]  uses  this  new  formulation  in  an  evaluation  of  the 
importance  of  shear  effects  in  scattering  from  the  sea  bed. 

1  FORMULATION 

Much  previous  work  on  acoustic  scattering  from  the  sea  bed  assumes  that  the 
sediment  can  be  treated  as  a  lossy  fluid,  in  which  case  scattering  at  the  in¬ 
terface  is  controlled  by  the  discontinuities  in  density  and  compressional  wave 
speed.  Volume  scattering  is  controlled  by  inhomogeneities  in  the  same  two 
parameters.  We  will  refer  to  these  controlling  factors  as  scattering  mecha¬ 
nisms.  Note  that  the  choice  of  compressional  wave  speed  as  a  bulk  parameter 
is  somewhat  arbitrary;  the  bulk  modulus  or  some  other  independent  parame¬ 
ter  could  be  used.  When  an  elastic  model  is  used  for  the  sediment,  additional 
degrees  of  freedom  appear  which  introduce  additional  scattering  mechanisms. 
In  particular,  one  must  now  consider  the  discontinuities  and  inhomogeneities 
in  shear  wave  speed  as  well  as  the  density  and  compressional  wave  speed. 
Furthermore,  the  scattering  mechanisms  may  generate  shear  waves  within 
the  sediment.  As  illustrated  by  Fig.  1,  this  leads  to  a  number  of  scattering 
mechanisms.  For  volume  scattering,  for  example,  one  mechanism  is  scatter- 
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ing  of  a  compressional  wave  by  a  density  fluctuation,  leading  to  conversion  to 
a  shear  wave.  This  shear  wave  can  be  converted  to  a  water-borne  compres¬ 
sional  wave  upon  encountering  the  interface.  If  the  shear  wave  is  polarized 
such  that  the  particle  motion  is  in  a  vertical  plane,  conversion  is  possible 
even  for  the  flat  interface.  As  a  result,  this  scattering  mechanism  is  evident 
in  a  flrst-order  perturbation  treatment.  Scattering  to  the  other  possible  po¬ 
larization  (with  particle  motion  in  a  horizontal  plane)  will  only  be  evident  at 
higher  order,  as  such  a  wave  suffers  total  internal  reflection  at  a  horizontal, 
flat  interface. 


1.1  Basic  Equations 

Figure  2  depicts  the  general  physical  situation  to  be  considered.  A  rough 
interface  described  by  the  random  function 

=  C(R)  (1) 

separates  a  lossless,  homogeneous  fluid  from  an  inhomogeneous  elastic  medium 
In  Eq.  1  and  following  equations,  upper-case  letters  will  be  used  to  denote 
two-dimensional  vectors  and  lower-case  letters  will  denote  three-dimensional 
vectors.  Thus,  the  three-dimensional  position  vector,  r,  is  decomposed  into 
its  transverse  and  vertical  components  as 


where 


r  =  (R,  0:3) 


(2) 


R=(xi,a:2)  •  (3) 

It  is  convenient  to  describe  wave  motion  in  terms  of  the  particle  displace¬ 
ment  field,  u(r).  This  field,  in  turn,  is  expressed  in  terms  of  scalar  and  vector 
potentials  [2] 


u(r)  =  V</i(r)  -I-  V  X  ^(r)  .  (4) 

These  potentials  describe  two  different  types  of  waves,  compressional  and 
shear,  supported  by  an  elastic  medium.  As  the  water  does  not  support  shear, 
the  vector  potential,  V’(r),  in  the  water  is  taken  to  be  zero. 

The  general  equation  obeyed  by  the  displacement  field  is 
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(5) 


(JpUi-{-di{W  •  u)+^9j  (/x(5-Uj+5jUi))  =  0  . 

i 

where  the  density,  /?,  and  the  Lame  parameters,  A  and  p,  can  be  functions 
of  position,  and  di  =  .  In  the  water,  p  is  taken  to  be  zero.  The  time 

dependence  of  all  field  quantities  is  assumed  to  be  exp(— iwt). 

The  boundary  conditions  on  the  interface  (xa  =  C(R))  ^.re  continuity  of 
the  normal  component  of  displacement 

u<^)-N  =  u-N  ,  (6) 

and  continuity  of  normal  tractions  for  the  stress  tensor,  crjj, 

O^ij  —  *  U  *h  ,  (*^) 


that  is, 

=  ,  (8) 

j  j 

where  the  superscript  ,  /,  denotes  the  fluid  medium  (water)  and  no  super¬ 
script  is  used  for  the  elastic  medium  (sea  bed  material).  The  vector,  N,  is 
the  (unnormalized)  interface  normal, 

N  =  -VC(R)  +  e3  ,  (9) 

where  63  is  the  unit  vector  in  the  vertical  direction. 

As  the  traction  boundary  condition,  Eq.  8,  is  actually  three  equation, 
there  are  a  total  of  four  equations  that  must  be  satisfied.  There  are  also 
four  unknowns:  the  scattered  compressional  field  in  the  water  and  the  three 
scattered  fields  in  the  sea  bed  material  (one  compressional  field  and  two  shear 
fields  having  horizontal  and  vertical  polarizations). 

1.2  Sea  Bed  Parameters 

As  mentioned  above,  other  parameters  can  be  used  instead  of  A  and  p.  In 
the  present  case,  compressional  and  shear  wave  speeds,  Cp  and  C(,  will  be 
employed  rather  than  the  Lame  parameters.  The  various  parameters  are 
related  as  follows; 
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A  =  p(cJ-2c^),  ii  =  pc].  (10) 

Subscripts,  p  and  t,  here  and  below  denote  compressional  (longitudinal)  and 
transverse  (shear)  waves,  respectively.  In  the  water,  Ct  and  p  are  taken  to  be 
zero. 

The  parameters  will  be  divided  into  two  classes:  those  representing  the 
average,  non-fluctuating  values  of  bulk  properties  such  as  density  and  wave 
speeds,  and  those  describing  the  random  fluctuation  of  bulk  properties  and 
interface  relief.  In  the  context  of  perturbation  theory,  the  zeroth-order  solu¬ 
tion  is  that  solution  corresponding  to  the  average  properties  with  a  planar 
interface.  The  first-order  correction  to  this  solution,  for  the  scattered  field, 
is  linear  in  the  fluctuating  parameters,  including  the  relief  of  the  interface. 

It  is  convenient  to  use  dimensionless  ratios  of  the  average  parameters.  For 
simplicity,  average  wave  speeds  and  density  are  simply  denoted  Cp,  C(,  and  p, 
although  these  symbols  were  earlier  used  for  the  position-dependent  quan¬ 
tities,  that  is,  the  average  value  plus  fluctuations.  This  simplified  notation 
will  be  employed  in  the  remainder  of  this  article.  Thus,  define 

a/ =  1  ,Gp  =  Cp/c/  ,at  =  ctlcf  ,ap  =  plpf  .  (11) 

The  trivial  ratio  aj  is  defined  for  convenience  in  later  expressions.  The  com¬ 
pressional  wave  speed  ratio,  Op,  is  the  ratio  of  the  average  sediment  com¬ 
pressional  wave  speed  and  the  water  compressional  wave  speed.  This  ratio  is 
complex,  with  imaginary  part  determining  absorption  loss.  The  shear  speed 
ratio,  at,  is  likewise  the  ratio  of  the  average  complex  sediment  shear  speed, 
and  the  compressional  wave  speed  in  the  water  and  the  density  ratio,  Up,  is 
the  ratio  of  the  average  sediment  mass  density  to  the  density  of  the  overlying 
water. 

Interface  roughness  will  be  described  by  thfe  spectral  density  function, 
$(>-)(K)  ,  where  K  is  the  two-dimensional  spectral  argument.  In  the  Born 
approximation,  this  is  a  sufficient  description  of  the  statistics  of  the  random 
interface.  In  other  approximations,  such  as  higher-order  perturbation  theory, 
additional  statistical  information  must  be  provided  by  specifying  higher-order 
statistics  or  by  assuming  Gaussian  statistics. 

A  description  of  the  statistics  of  volume  inhomogeneities  is  given  by  a  ma¬ 
trix  of  spatial  (3-dimensional)  cross-spectra.  These  are  Fourier  transforms 
of  cross-correlation  functions  for  the  relative  fluctuations  Cp  =  Apfp,Cp  = 
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Acpfcp,et  =  Actfct.  These  relative  fluctuations  will  be  collectively  desig¬ 
nated  tjs  with  P  =  p,p^t.  In  this  notation,  the  cross-spectra  are  $|gjg/(K,(2) 
where  K  and  q  are  the  horizontal  and  vertical  components  of  the  three- 
dimensional  spectral  argument.  Again,  these  spectra  provide  a  sufiicient 
description  of  the  inhomogeneities  in  the  Born  approximation. 

2  THEORETICAL  RESULTS 

The  end  result  of  the  formulas  to  be  presented  is  the  scattering  coefficient, 
(7,  or  scattering  cross  section  per  unit  area  per  unit  solid  angle.  The  term 
“cross  section”  will  be  used  for  brevity,  even  though  this  quantity  is  dimen¬ 
sionless.  The  quantity  lOlogcr  is  commonly  called  the  scattering  strength 
and  is  expressed  in  decibels.  The  scattering  cross  section  characterizes  the 
frequency-angular  distribution  of  the  mean  intensity  of  field  fluctuations  in 
the  Fraunhofer  zone  (or  far  zone,  relative  to  a  scattering  surface  or  volume). 
In  the  far  zone,  the  scattered  acoustic  wave  can  be  approximated  as  a  plane 
wave.  Likewise,  it  is  implicit  in  the  definition  of  scattering  strength  that  the 
scattering  interface  is  in  the  far  zone  of  the  source,  so  that  the  incident  field 
can  also  be  approximated  by  a  plane  wave.  The  scattering  cross  section  is 
assumed  here  to  be  the  sum  of  contributions  from  interface  roughness  and 
volume  inhomogeneity.  In  the  first-order  approximation  used  here,  this  will 
be  the  case  if  the  interface  relief  is  uncorrelated  with  the  volume  inhomo¬ 
geneities  [3,  10,  14,  20]. 

cr(K„K.)  =  cr,(K„K..)  -f  a„(K„K,)  (12) 

Expressions  for  these  two  components  are  given  without  derivation  in  the 
two  following  subsections.  The  appendices  give  outlines  of  the  derivation  of 
(Tr(Ki,K,)  and  cri,(Ks,Kj)  using  unified  notation.  In  addition,  the  volume 
scattering  formalism  of  Ref.  [13]  is  generalized  to  include  a  wider  class  of 
scattering  mechanisms.  The  formulas  presented  in  this  section  are  written  in 
a  form  convenient  for  coding  in  higher-level  matrix-oriented  computational 
language. 
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2.1  Relevant  Plane  Waves 

It  is  convenient  to  present  a  solution  of  the  first-order  scattering  problem 
in  terms  of  solutions  of  the  unperturbed  wave  equations,  that  is,  the  plane 
waves  propagating  in  the  up  and  down  directions, 

exp(ik*  •  r) ,  a  =  /,  p,  t  , 

with  the  subscripts  a  denoting  the  wave  vectors  for  compressional  waves  in 
the  water,  compressional  waves  in  the  sea  bed,  and  shear  waves  in  the  sea 
bed,  respectively.  The  superscripts  +  and  —  denote  the  direction  of  the  wave 
propagation,  up  and  down,  respectively. 

It  is  also  appropriate  to  define  the  general  three-dimensional  wave  vec¬ 
tors  through  the  transverse  components  of  the  incoming  and  outgoing  wave 
vectors 

kf(K)  =  (K,±i,./.(K)),a  =  /,  p,  i  (13) 

The  unperturbed  wave  equations  obeyed  by  the  potentials  determine  the 
relation  between  the  transverse  and  vertical  components  of  the  wavevectors: 

i/a(K)  =  ^a-2  -  lOjk)  ,a  =  fyp,t  .  (14) 

The  speed  ratios,  Cq,  were  defined  earlier. 

The  following  unit  vectors  specifying  the  directions  of  propagation  are 
also  useful. 

e±(K)  =  kJ(K)/J:.  ,  (15) 

where 

u 

kc  =  —  ,Q  =  /,p,t,  (16) 

Cq 

are  the  wave  numbers,  respectively,  of  compressional  waves  in  water,  com¬ 
pressional  waves  in  the  sea  bed,  and  shear  waves  in  the  sea  bed. 

There  are  two  plane-wave  shear  polarizations  to  consider,  both  having 
particle  displacement  normal  to  the  direction  of  propagation.  The  shear 
wave  having  particle  displacement  in  the  direction 

e±(K)  =  e3xK/A^  =  (-A'2/AMih/A;0),  (17) 

will  be  referred  to  as  “horizontally  polarized”  as  the  particle  displacement  is 
in  a  horizontal  plane.  In  this  case,  the  polarization  vector  does  not  depend 
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upon  whether  propagation  is  upward  or  downward,  but  the  superscripts  ± 
are  used  for  consistency  with  the  vertical  polarization  case.  A  plane  shear 
wave  having  transverse  wave  vector,  K,  and  particle  displacement  in  the 
direction  of  the  unit  vector 

e±(K)  =  kf  (K)  X  e±(K)/fct  =  TKatUt{K)IK  +  e^K/h  ,  (18) 

will  be  referred  to  as  “vertically  polarized”,  as  the  displacement  is  in  a  vertical 
plane.  In  this  case,  the  polarization  vector  depends  upon  whether  the  wave 
is  propagating  upward  or  downward.  Scattering  by  interface  roughness  and 
inhomogeneities  within  the  seabed  causes  conversion  between  the  wave  types 
defined  above.  For  example,  while  incidence  of  a  plane  compressional  wave 
in  the  water  on  a  flat  interface  only  gives  rise  to  a  vertically  polarized  shear 
wave,  scattering  by  interface  roughness  causes  conversion  of  compressional 
energy  in  the  water  to  both  the  vertical  and  horizontal  shear  polarizations. 
Scattering  by  inhomogeneities  within  the  seabed  causes  conversion  between 
all  three  wave  types. 

The  above  expressions  for  wave  vectors  and  polarization  vectors  are  writ¬ 
ten  in  terms  of  a  general  transverse  wave  vector,  K,  for  later  use  in  the 
theoretical  development.  In  computing  the  scattering  cross  section,  however, 
one  must  take  K  =  Ki,^.  In  this  case,  it  is  convenient  to  define  the  transverse 
components  in  terms  of  the  angular  coordinates  of  interest: 

K(0,  </»)  =(A:/ cos  0  cos  <^,  cos  0  sin</))  ,  (19) 

and 

K,>  =  K(^,,„0,-,),  (20) 

where  Oi  and  0s,  are  the  grazing  angles  of  the  incident  and  scattered  acoustic 

waves,  respectively,  and  4>s  and  4>i  are  the  azimuthal  angles  of  the  incident 
and  scattered  waves,  respectively.  The  unit  vectors  can  then  be  expressed  as 
follows: 

e^(K)  =  a  a[cos0  COS(l>,COs6  Sm(t>,±Ua{K)]  ,  (21) 

ef(K)  =  [-sin^i!),cos(jf>,0]  ,  (22) 

e*(K)  =  at[T^'t(K)cos^,  =Fr't(K)sind>,  COS0]  ,  (23) 

with 

i/o(K)  =  —  cos^  9.  (24) 
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In  Eqs.  21-24,  subscripts  i  and  s  must  be  added  to  K  and  to  the  angles 
,  0  and  (f>  as  appropriate. 


2.2  Scattering  due  to  Interface  Roughness 

The  derivation  of  the  results  presented  here  is  given  in  Appendix  B.  The 
scattering  cross  section  for  interface  roughness  is  of  the  form 

(7,(K„K,.)  =  |///(K„K,.)P$<’-HK.-Ki)  ,  (25) 


where  $^'"1  is  the  roughness  spectrum  defined  in  Eq.  50,  and  .^/(Ks,K,)  is 
the  first  element  of  the  column  matrix 


//(K„K,.) 


Hj{Ks,Ki)  \ 
^p(K„K,) 

HhiK^Ki)  / 


(26) 


As  explained  in  Appendix  B,  the  elements  of  Hf{Ks,  Kj)  determine  the  plane- 
wave  spectra  of  energy  scattered  into  the  water,  into  compressional  waves  in 
the  sea  bed,  and  into  vertically  and  horizontally  polarized  shear  waves  in  the 
sea  bed.  This  matrix,  in  turn,  is  computed  as  follows: 

H{K,,Ki)  =  A:Jyi(K,)[p(3)(K,)]-iB(K„KJZ)o(K.)  .  (27) 

In  this  equation,  Do(K,)  is  a  five-row  column  vector  comprised  of  the  four 
transformation  (reflection-transmission)  coefficients  contained  in  the  column 
vector,  S'o(K,),  and  supplemented  with  unity  in  the  last  row. 


Do{Ki) 


So{Ki)  \ 
\  1  / 


(28) 


An  expression  for  computation  of  5'o(K,)  will  be  given  later.  Matrix  B(Ks,  K-) 
is 


B(K„K,)  =  kf{K,,-K„)E^^\Ki)+kj\Ks2-Ki2)E^^\Ki)-Y,{Ki)E^^\Ki) . 

(29) 
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Note  that  the  factors  kJ^{Ksn  —  Kin),  n  =  1,2,  are  dimensionless  and  can 
be  expressed  conveniently  in  terms  of  the  angles  defined  in  Eq.  19.  The 
matrices  Yi  and  V2  are 


^i(K) 


Uf{K)  0  0  0  \ 

0  i/p(K)  0  0 

0  0  i/t(K)  0 

0  0  0  ut{K)  / 


(30) 


/  uj{K)  0  0  0  0  \ 

0  -i^p(K)  0  0  0 

y2(K)  =0  0  -utiK)  0.0  .  (31) 

0  0  0  -ut{K)  0 

V  0  0  0  0  -ujiK) 

There  are  three  matrices  E^"1(K),  defined  for  each  of  the  coordinate  indices, 
n  =  1,2,3 

/  Pjr’(K)  P,‘,”’(K)  Fir’(K)  p£'(K)  q1”>(K)  \ 

F(“)rKi=  I  m', 

'  P)3"'(K)  P^”>(K)  Pi3">(K)  pI^\K}  I  QP(K)  ’ 

\  pj:>(K)  p<r>(K)  p<;>(K)  pi;'(K)  q1"'(k)  j 

where  the  dashed  vertical  line  in  Eq.32  separates  it  into  the  the  4x4  matrices, 
p(”)(K),  and  the  column  matrices,  (5^”^(K).  Note  that  P(^)(Ks)  appears  in 


Eq.  27.  The  elements  of  P("^(K)  and  (5*"^(K)  are 

P)”>(K)  =  i„„,  m  =  1,2,3  ,  (33) 

Pj;>(K)  =  e+„(K)  ,  (34) 

^pTnC^)  ”  ^^["“^77171  +  2(at/ap)  {Synn  “  ^pm^pn)]  7  HI  =  1,2,3  ,  (35) 

PM(K)  =  -aj‘e;„(K)  ,  (36) 

P£i)».(K)  = -Ojlermer.,).),  +  er,er.,*)J  •  "!=  1,2,3  ,  (37) 

P£i)4(K)  =  ,  (38) 
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(3<;>(K)  =  ,  m  =  1,2,3  ,  (39) 

Oi”’(K)  =  e7„(K)  ,  (40) 

where  5nm  is  the  Kronecker  delta  function  and  the  unit  vectors  defined  earlier 
are  enaployed. 

The  transformation  matrix,  5o(Kj),  is  computed  as  follows: 

5o(K0  =  .  (41) 

The  four  elements  of  this  matrix  are  denoted 


5o(K.) 


W,(K,)  \ 
WAK,) 
W,(K,)  ’ 

Wh(KA  ) 


and  are  the  in- water  reflection  coefficient  for  compressional  waves,  and  the 
transmission  coefficients  for  compressional  waves  and  the  two  shear  w'aves  in 
the  sea  bed  (see  Appendix  B).  Alternative  expressions  for  these  coefficients 
are  given  in  Ref.  [2].  Note  that  Wh{Ki)  =  0.  These  coefficients  are  also 
required  for  the  volume  scattering  computation. 


2.3  Scattering  due  to  Volume  Inhomogeneity 

The  scattering  cross  section  for  volume  inhomogeneity  is  of  the  form 

(43) 

where  Im  designates  the  imaginary  part  of  the  complex  quantity  and  $^73/ 
is  the  is  the  matrix  of  cross  spectra  for  volume  inhomogeneities  defined  in 
Eq.  51.  Thus,  the  (3  and  (3'  sums  run  over  the  three  types  of  inhomogeneity: 
density,  compressional  wave  speed,  and  shear  wave  speed  (/?,/?'  =  p,p,t). 
Generally,  these  fluctuations  are  expected  to  be  strongly  correlated,  thus  the 
^  ^  (3'  cross-terms  are  important.  The  t]  and  rj'  sums  run  over  the  four 
types  of  wave  conversion  caused  by  volume  scattering:  1.  compressional  to 
compressional  (77  =  pp),  2.  shear  to  compressional  (77  =  pt),  3.  compressional 
to  shear  (77  =  tp),  and  4.  shear  to  shear  (77  =  tt).  Thus,  we  can  put 


<7‘t/(Ks,  K,)  — 


((qu  +  q;>)/2) 

(9173  - 
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T]  =  aa', 

where  a  and  a'  run  over  the  two  types  of  waves  (a,  a'  =  p,  t).  As  noted  ear¬ 
lier,  volume  inhomogeneities  cause  conversion  to  both  shear  polarizations, 
but  upgoing  shear  waves  with  horizontal  particle  motion  suffer  total  reflec¬ 
tion  at  the  interface  and  do  not  contribute  to  the  scattering  cross  section  in 
first  order.  The  scattering  vectors  for  the  four  relevant  types  of  conversion 
corresponding  to  different  channels  of  scattering  are 

q„„.  =  k+(K,)  -  k;,(K,)  (44) 

These  vectors  give  the  change  in  wave  vector  for  the  corresponding  conversion 
and  appear  as  arguments  in  the  cross  spectra  in  Eq.  43.  All  of  these  vectors 
have  the  same  transverse  components 

q„  =  (K^-Ki,q„3)  .  (45) 

The  vertical  components  are 

The  coefficients  are 

=  W-qD-riP  (4"^) 

where 

Waa'  ^  Wa(K5)W«^(Ki),  a,  a'  =  p,t  (48) 

with  Wi  =  Wy.  The  coefficients  Dap  are  elements  of  the  3-column  matrix 

D  =  {D,\D^\D,) 

where 

-2a]a-'^  +  2a]hl^- 

^pv  4"  ^pv^pt 
byp  hyphpf 

byy  a^b 

VV  btt  T  dtbyibt 

D,  =  {2,0,0, Of 


(46) 
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Dt  =  — 2aj 


Here 


Ka'  =  K{K,)-k-,{Ki)/k}  (49) 

where  a,  a'  =  p,t,v  and  the  wave  vectors  kJ(K)  are  defined  in  Eq.  13. 


3  CALCULATIONS  AND  DISCUSSION 


3.1  Model  Parameters 

To  illustrate  some  of  the  effects  of  shear  elasticity  on  sea  bed  scattering, 
calculations  will  be  presented  for  a  single  set  of  parameters,  with  average 
physical  parameters  appropriate  to  sedimentary  rock  [5].  These  are  density 
ratio,  Gp  =  2.5,  compressional  speed  ratio,  Op  =  2.3  —  0.004i,  and  shear 
speed  ratio,  at  =  1.3  —  0.22i.  The  water  sound  speed  is  taken  to  be  1500 
m/s  and  the  frequency  is  taken  to  be  1  kHz.  As  noted  later,  the  roughness 
and  volume  fiuctuation  parameters  used  are  such  as  to  yield  essentially  no 
frequency  dependence. 

The  roughness  spectrum  is  taken  to  be  of  the  form 


(2/7r)A(’') 

{hi  + 


(50) 


Equation  50  assumes  that  roughness  statistics  are  stationary  and  isotropic, 
with  a  spectrum  that  obeys  a  power  law  for  Kq  .  The  parameter 

Kq  is  an  inverse  correlation  scale  for  roughness  in  the  transverse  directions. 
The  roughness  spectral  parameters  were  assigned  the  following  values:  Kq  = 
10~^m“^  and  =  10~^.  The  form  of  the  spectrum  and  parameter  Kq 
were  chosen  so  that  the  scattering  cross  section  is  essentially  independent 
of  frequency  except  for  scattering  very  close  to  the  specular  direction.  The 
parameter  is  dimensionless  and  was  chosen  arbitrarily  to  yield  scattering 
strengths  that  are,  except  near  the  specular  direction,  within  the  likely  region 
of  validity  of  the  Born  approximation.  The  factor  2/7r  is  included  in  Eq.50 
for  convenience  [17]. 
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The  volume  inhomogeneity  spectra  were  taken  to  be  of  the  form 


(K,9) 


(A'2 +  52  +  52)3/2- 


(51) 


with  go  =  10"^m“^,  =  10~®.  These  spectra  correspond  to  isotropic 

uncorrelated  inhomogeneities.  The  form  of  the  spectra  and  parameter  go 
are  chosen  so  that  the  scattering  cross  section  is  essentially  independent  of 
frequency  for  all  scattering  directions.  The  parameter  is  dimensionless 
and  was  chosen  under  the  same  rationale  as  A^^^  with  a  factor  1  / (27r)  included 
for  convenience  [17]. 

Because  of  the  Kroneker  delta  used  in  Eq.  51,  all  the  nondiagonal  el¬ 
ements  of  the  correlation  matrix  vanish.  In  this  “uncorrelated”  case,  the 
volume  component  of  the  scattering  cross  section,  according  to  Eq.43,  will 
be  simply  the  sum  of  three  components 


=  (52) 

where  cr^j  correspond  to  contributions  of  fluctuations  of  the  three  different 
bulk  parameters  (/?  =  p,Cp,Ct).  Calculations  are  presented  below  for  each  of 
the  terms  in  Eq.  52.  For  example,  Up  is  considered  for  the  case  in  which  there 
are  density  fluctuations  but  no  compressional  or  shear  speed  fluctuations. 


3.2  Numerical  Results 

The  parameters  specified  above  were  used  in  computing  angular  dependencies 
of  the  backscattering  strength.  The  results  of  these  calculations  are  shown 
in  Fig.  3  which  shows  the  relative  contributions  of  different  types  of  per¬ 
turbation  for  the  elastic  case  and  compares  them  to  the  fluid  case.  As  three 
different  types  of  volume  fluctuation  are  considered,  this  comparison  involves 
a  set  of  seven  different  curves  that  includes  three  pairs  comparing  the  elastic 
and  fluid  cases  for  scattering  by  roughness,  by  density  fluctuations,  and  by 
compressional  wave  speed  fluctuations.  The  seventh  curve  is  for  scattering 
by  shear  wave  speed  fluctuations,  which  have  no  fluid  counterpart. 

The  results  show  that  for  sedimentary  rock  the  effects  of  shear  elasticity 
on  both  roughness  and  volume  scattering  are  very  strong.  As  the  average 
parameters  are  taken  from  Essen  [5],  the  roughness  scattering  results  are 
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equivalent  to  his.  Figure  3  shows  that  roughness  backscattering  strength 
exhibits  marked  dips  at  grazing  angles  immediately  below  the  critical  angles 
for  shear  and  compressional  waves  (39.7°  and  64.2°,  respectively).  The  inclu¬ 
sion  of  lossy  elasticity  is  seen  to  decrease  the  roughness  scattering  strength 
for  certain  angles  compared  to  the  fluid  case  in  which  the  shear  speed  is  set 
to  zero.  But,  at  the  same  time,  it  somewhat  enhances  scattering  for  other 
angles  (in  the  40°  -  60°  range).  Clearly,  treating  sedimentary  rock  as  a  fluid, 
albeit  a  massive  and  stiff  fluid,  is  a  poor  approximation. 

Figure  3  shows  that,  in  contrast  to  the  case  of  roughness  scattering, 
elasticity  greatly  enhances  volume  backscattering  in  rock  for  grazing  angles 
smaller  than  the  shear  wave  critical  angle  for  all  the  types  of  volume  fluctua¬ 
tions,  but  can  cause  a  great  reduction  in  scattering  from  compressional  speed 
fluctuations  for  grazing  angles  in  the  40°  -  60°  range.  Elasticity  effects  for 
volume  scattering  are  much  larger  than  the  effects  for  roughness  scattering. 

As  these  results  show,  the  effects  of  elasticity  on  roughness  and  volume 
backscattering  for  rock  are  strong  and  complicated.  A  more  detailed  numer¬ 
ical  study  of  the  effects  of  elasticity  for  different  sea  bed  types,  different  sets 
of  statistical  parameters,  and  different  types  of  geometry  (both  monostatic 
and  bistatic  cases)  will  be  given  in  a  subsequent  article  [17]. 

4  CONCLUSIONS 

First-order  perturbation  theory  has  been  used  to  develop  expressions  for  the 
scattering  strength  of  elastic  sea  beds.  Both  roughness  and  volume  compo¬ 
nents  of  scattering  are  included.  The  roughness  component  of  the  scattering 
strength  is  related  to  the  spatial  two-dimensional  spectrum  of  roughness.  The 
volume  component  is  expressed  through  three-dimensional  cross-spectra  for 
spatial  fluctuations  of  density,  compressional  and  shear  wave  speed.  Shear 
effects  on  both  roughness  and  volume  scattering  are  complicated  and  strong, 
at  least  for  the  case  of  backscattering  from  the  sedimentary  rock  taken  as  an 
example  for  calculations.  Generally,  an  assessment  of  shear  effects  for  various 
sea  bed  types  requires  more  detailed  numerical  study,  using  available  data 
on  sea  bed  parameters  including  their  statistics. 
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Appendix  A:  Scattering  Amplitudes 

Regardless  of  the  complexity  of  the  random  sea  bed,  whether  scattering  is 
due  to  roughness  or  volume  inhomogeneity,  the  scattered  field  in  the  water 
above  the  highest  point  on  the  interface  can  be  expressed  as  a  superposition 
of  plane  waves.  In  the  Born  approximation,  plane-wave  superpositions  are 
also  used  to  describe  the  fields  in  the  sediment.  If  the  incident  field  is  a  plane 
wave  of  unit  amplitude  having  wave  vector  with  transverse  component,  K,, 
then  the  scattered  potential  in  the  water  is 

lO  T  ^ 

MK,K,)  exp(ik+(K)  ■  r)  ,  (53) 

where  A/(K,K,)  is  the  scattering  amplitude  [13],  related  to  the  transition 
matrix  (“T  matrix”)  [33,  29]  as  follows: 

yl,(K,K,)  =  fc/i./(K)r(K,Ki)  .  (54) 

The  potential,  ^/(r),  does  not  include  the  incident  field  but  does  include 
all  other  portions  of  the  field  (coherent  and  incoherent).  The  scattering 
amplitude  obeys  the  reciprocity  relation 

Af{K,K,)  =  Af{-Ki,-K)  .  (55) 

Similarly,  the  total  scalar  potential  in  the  sea  bed  is  expressed  as 

^  /  kjUp{K)  exp(ik;(K)  •  r)  ,  (56) 

and  the  vector  potential  in  the  sediment  is 

j2  t  ^ 

/  kjUt{K) 

where 

A,(K,K,)  =  e^(K)A4K,Ki)  -  e;(K)/l,(K,Ki)  (58) 
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When  the  in-water  scattering  amplitude  has  been  found,  all  the  basic  charac¬ 
teristics  of  the  scattered  field  can  be  determined.  For  example,  the  coherent 
reflection  coefficient  is  defined  by  the  following  equation  for  the  averaged 
scattering  amplitude,  or  its  coherent  component 

<  Af{K,  Ki)  >=  5{K  -  Ki)Wf{Ki)i^f{Ki).  (59) 

Generally,  the  coherent  reflection  coefficient,  W /(K,),  is  not  equal  to  the 
zeroth-order  reflection  coefficient,  but  this  distinction  is  not  significant  in 
the  present  problem. 

The  scattering  amplitude  will  be  used  to  obtain  the  scattering  cross  sec¬ 
tion  through  the  equation  [33,  31,  29,  1] 

<  A' (K„  k,)a';(k;,  k,)  >=  5(k,-k' )<7(k„  k,)  ,  (eo) 

where  A'^(Ks,KJ  =  Af{Ks,Ki)—  <  A/(Ks,K,)  >  is  the  incoherent  com: 
ponent  of  the  scattering  amplitude  and  it  is  assumed  that  the  statistics  of 
the  random  sea  bed  are  stationary  in  the  the  two  transverse  (horizontal) 
coordinates.  The  water  column  is  assumed  to  be  homogeneous,  supporting 
plane  waves,  and  the  scattering  cross  section  is  assumed  to  be  the  sum  of 
contributions  from  interface  roughness  and  volume  inhomogeneity. 

Appendix  B:  Interface  Scattering 

The  incoherent  roughness  scattering  amplitudes  in  the  Born  approximation 
have  the  general  form 

<(K„K,)  =  iif„(K„K,)Z(K,-K,)  ,a  =  /,p,u,fi  ,  (61) 

where  Z(K)  is  the  Fourier  transform  of  the  interface  relief  function,  Eq.  1. 
Then  Eq.  60  leads  directly  to  the  equation  for  the  roughness  scattering  cross 
section,  Eq.  25.  Thus,  the  purpose  of  this  appendix  is  to  outline  the  deriva¬ 
tion  of  the  roughness  scattering  amplitude  and,  particularly,  Eq.  27.  As 
discussions  of  this  problem  have  been  given  by  others  [4,  5],  this  appendix 
will  only  give  a  brief  outline  that  emphasizes  the  physics  and  establishes  no¬ 
tation  convenient  for  computation  and  consistent  with  the  volume  scattering 
case. 
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The  scattering  amplitude  is  determined  by  imposing  the  boundary  condi¬ 
tions,  Eqs.  6  and  8,  on  the  fields.  The  derivative  operations  that  yield  the  dis¬ 
placement  and  shear  are  readily  implemented  in  the  plane- wave  expansions, 
where  derivatives  become  multiplication.  For  example,  the  displacement  in 
the  water  is 

=  i  /  ^  e+(K)A;(K,  K,)  exp(ik+(K)  ■  r)  ,  (62) 

and  the  displacement  due  to  shear  in  the  sea  bed  is 

jO 

uW(r)  =  ia,-‘/ ^  |e;(K)A„(K,K.)  +  e;(K)Ai(K,K,)lexp(ikr(K)T) 

(63) 

The  additional  derivatives  needed  to  form  the  stress  tensor  introduce  ad¬ 
ditional  wave  vector  factors,  which  are  taken  as  unit  vectors  multiplied  by 
wave  numbers.  The  derivative  factors  comprise  the  matrix  given  in  Eq.  32. 
The  first  column  of  pertains  to  the  in-water  scattering  amplitude, 

while  the  second,  third,  and  fourth  pertain  to  the  p,  u,  and  h  fields  in  the  sea 
bed  material.  The  fifth  column  pertains  to  the  incident  field.  The  boundary 
conditions  require  that  the  differences  in  normal  displacement  and  normal 
tractions  vanish  across  the  boundary,  accordingly,  the  matrix  con¬ 

tains  appropriate  sign  reversals.  The  first  three  rows  (m  =  1,2,3)  of  £'("')(K) 
express  the  derivatives  needed  for  the  stress  tensor,  amn,  and  the  last  row  ex¬ 
presses  the  derivatives  needed  for  the  nth  component  of  displacement.  That 
is,  considering  the  plane-wave  expansion 

=  JcPK  E^^\K)Yf\K)U  ,  (64) 

with 

(  A/(K,KJ  exp(ik^(K)  •  r)  ^ 
x4p(K,K,)  exp(jk“(K)  •  r) 
A,(K,KJexp(ikr(K)T)  ,  (65) 

A/i(K,  K,)  exp(ik7  (K)  •  r) 

V  kjUj{K)5{K  -  Ki)  exp(iko  (K)  •  r)  / 

the  mth  row  of  =  1,2, 3)  is  proportional  to  the  water-sea  bed  differ¬ 

ence  in  the  mn-component  of  stress,  and  the  4th  row  is  proportional  to  the 
water-sea  bed  difference  in  the  nth  component  of  the  displacement  vector. 
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The  boundary  conditions  can  now  be  written  as 


+  =  0  ,  (66) 

with  X3  =  C(K)-  To  obtain  the  zeroth-  and  first-order  solutions,  the  scattering 
amplitudes  are  expressed  as  a  sum  of  zeroth-  and  first-order  terms,  the  former 
being  the  flat-surface  solution  and  the  latter  being  linear  functionals  of  C(R)- 

A.(K,K,)  =  /1!,“>(K,K,)  +  <(K,K,.)  .  (67) 

Also,  the  exponentials  used  are  expanded  to  first  order  in  C(R)  with  X3  = 

C(R) , 

exp(ikJ(K)  •  r)  =  exp(iK  •  R)[l  ±  iA:/i/a(K)C(R)  -f  •  •  •]  .  (68) 

Finally,  a  Fourier  transform  of  Eq.  66  with  respect  to  R  is  taken  and  the 
definition 

Z(K)  =  (2^  /  d^Rexp{-iK  ■  R)C(R)  (69) 

is  used,  and  the  derivatives  of  C(R)  Eq.  66  are  treated  using 

z'KZ(K)  =  ^  J  d^Rexvi-iK  ■  R)VC(R)  .  (70) 

The  zeroth-order  terms  yield 

P(^)(K,)A(°)(K„  K,)  -f  Q(")(K0^(K,  -  K.)  =  0  ,  (71) 

where  A^°)(Ks,K,)  is  a  column  matrix  comprised  of  the  zeroth-order  scat¬ 
tering  amplitudes.  Equation  71  is  equivalent  to  Eq.  41. 

The  equation  for  first-order  terms  is 

ikpiK,,Ki)Do{K,)Z{K,  -  K,)  -  F(^)(K,)yriA'(K„K,)  =  0,  (72) 

where  P(Ks,Ki)  is  given  by  Eq.  29  and  A'(Ks,Kj)  is  a  column  matrix 
comprised  of  the  first-order  approximation  to  the  incoherent  scattering  am¬ 
plitudes.  Equation  72  is  equivalent  to  Eqs:  27  and  61. 
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Appendix  C:  Volume  Scattering 

This  appendix  outlines  first-order  perturbation  theory  (Born  approximation) 
for  scattering  from  an  inhomogeneous  elastic  half  space  with  planar  boundary 
and  fluid  above.  This  treatment  is  the  same  in  principle,  but  less  detailed 
than  that  given  by  Ivakin  [13].  This  appendix  uses  different  notation  than 
Ref.  [13],  and  a  more  general  expression  for  the  bistatic  cross  section  is 
developed  here,  in  a  form  convenient  for  computation. 

Consider  small  fluctuations  of  elastic  medium  parameters  and  correspond¬ 
ing  corrections  to  the  displacement  field  and  potentials,  making  correspond¬ 
ing  change  in  the  wave  equations 

^  p  +  p',X  +  +  +  +  + 

This  gives,  for  the  zero-order  potentials,  the  following  well-known  wave  equa¬ 
tions  for  a  homogeneous  elastic  medium 

(V^  +  kl)4>  =  0. 

(V" + = 0, 

and  the  following  equation  for  the  first-order  field  perturbation  in  the  dis¬ 
placement 


pJ^u'  /!>CpV(V  •  u')  —  pcj  V  X  V  X  u'  =  f , 
where  f  is  a  vector  with  components 


fi  =  -p'u'^ui  -  d,  (A'(V  •  u))  -  ^  dj  {p'idiUj  +  djUi)) ,  i,  j  =  1,2, 3.  (73) 

j 

This  corresponds  to  the  following  pair  of  first-order  equations  for  the  per¬ 
turbed  potentials: 


p{u^  -F  cJV2)VV'  =  V  •  f 
p{lj^  +  cy^)V  xVx^'  =  Vxf. 
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A  solution  of  these  equations  is 

= -^  /  Gp(r  -  r')V  •  f(r'Kr', 
lAW  /  G,(r  -  r')V  X  f(rVr', 

where  Gc,{r){a  =  p,t  )  are  free-space  Green’s  functions,  which  can  be  ex¬ 
pressed  as  expansions  in  plane  waves  propagating  up  toward  the  interface 

G«(r)  =  -  exp(ii.r)  =  -^  /  exp  (ik:(K)  •  r)  . 

Then  the  scattered  fields  which  are  incident  on  the  interface  can  be  also 
expressed  in  terms  of  plane  wave  expansions 

J  kfiyp{K)  •  r)  , 

~  /  kjUt{K)  ■  J*)  > 

where 


Z7T 


4  =  r:j^kJv.f), 


puj 

in 

pLU 


tn 

a:  =  --Tfk,(V  X  f), 


■^k(Q)  =  (27r)  ^  j  Q{r)exp{-ik- r)(fr.  (74) 

Transmission  of  these  waves  through  the  flat  interface  is  governed  by  the 
expression 


_  a;(k.k,) 


+  W+(K) 


e,(K).A;(K,K,) 

r,{K) 


5 
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where  the  are  the  transmission  coefficients  for  the  upward  direction  and 
are  related  to  downward  direction  coefficients,  Wa,  as  follows 

^a/^o,  =  apWajuj,  a  =  p,t. 

Note  that  Wt  =  Wy. 

Then,  using  the  zero-order  solution  in  Eq.  73,  one  obtains  the  following 
expression  for  the  scattering  amplitude 


^/(K,K,)  =  (75) 

with  q„  as  defined  by  Eqs.44-46  and  dr,/3  as  defined  by  Eqs.  47-49.  Using  Eqs. 
74,75,  the  expression  for  the  cross  section,  Eq.  43,  can  be  readily  obtained. 
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